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ABSTRACT 

We compared the total mass density profiles of three different types of galaxies using 
weak gravitational lensing: (i) 29 galaxies that host quasars at z ~ 0.32 that are in a 
post-starburst (PSQ) phase with high star formation indicating recent merger activity, 
(ii) 22 large elliptical galaxies from the SLAGS sample that do not host a quasar at 
z ~ 0.23, and (hi) 17 galaxies that host moderately luminous quasars at z ^ 0.36 
powered by disk instabilities, but with no intense star formation. On an initial test 
we found no evidence for a connection between the merger state of a galaxy and the 
profile of the halo, wit h the PSQ profile compa rable to that of the other two samples 

(20141 study of moderately luminous quasars 
the two quasar samples, we combined these 


and consistent with the Leauthaud et al. 
in COSMOS. Given the compatibility o: 


and found no evidence for any connection between black hole activity and the dark 
matter halo. All three mass profiles remained compatible with isothermality given the 
present data. 
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1 INTRODUCTION 

The mass of a black hole (BH) has been shown to tightly 
correlate with the size of the stellar bulge of its host galaxy 
and even tighter with its stellar velocity dispersion 


[rarese &i Merritt 2000 Gebhardt et al.||2000 l, implying an 
inextricable link between galaxy evolution and BH forma¬ 
tion ( [Englmaier fc Shlosman|2000[ [Emsellem et al.||2014[ ) . 

Efforts to construct a single unifying theory that can 
explain the formation and fueling of BHs were set back with 


the advent of the Hubble Space Telescope (HST) (e.g., Bah- 
|call et ah 1997 Disney et al.||19S)5 K With the improved an¬ 
gular resolution of a space-based telescope and with deep 
spectroscopy, BHs were found indiscriminately in old, early- 
type elliptical galaxies and younger, smaller, spiral galaxies 
(e.g. Letawe et al.|2008 2007 Jahnke et al.|2007 1. As large 
elliptical galaxies with brighter central quasars are more 
common at higher redshifts than spiral galaxies with less 
luminous BHs, and vice versa at lower redshifts, it was sug¬ 
gested that the average cosmic quasar activity for a given 
galaxy dropped as the Universe evolved to lower redshift 


(e.g Dunlop et al.|2003 McLure et al.|1999 Heckman et al. 

|2004D . 

This “Cosmic Downsizing” suggested that there existed 
multiple fueling mechanisms that could sustain BHs, lead¬ 


ing to a two-mechanism paradigm in BH formation (Barger 


|et al.|2005||Hasinger et al.|2005p . The first mechanism to pro¬ 
duce a BH occurs earlier in the Universe where the merger of 
two galaxies can instigate an intense burst of star formation 
and strong inflows, triggering the ignition of a quasar. The 
quasar phenomenon, the luminous counterpart of a super- 
massive BH, is an important stage of a galaxy’s life, since 
the resulting feedback blows out any local dust, quenching 
star formation and ultimately evolving the galaxy onto the 


red sequence (|Di Matteo et al. 2005 Springel et al. 2005 
[Hopkins et al.pOOOp . The second mechanism occurs in spi¬ 
ral galaxies at lower redshifts where a large central bar can 
contain sufficient gas to continually accrete and fuel a BH 
i jHopkins fc Hernquist|2006||2009| ). This mechanism results, 
however, in much less luminous quasars. 

In this letter we compared the effect these two mecha¬ 
nisms had on the density profile of quasar host galaxies using 
weak gravitational lensing. To study the former mechanism, 
i.e. that BHs are triggered by merging events, we looked at 
galaxies that showed recent signs of intense star formation, 
as it is thought that any major merger will instigate the for¬ 
mation of large quantities of stars. The HST snapshot sur¬ 
vey of 29 Post-Starburst Quasar (PSQ) showed clear signs 


of recent intense star formation (Gales et al. 2011 20131 


and is was suited to this study. This point, during or just 
after a merger could theoretically have produced a steeper 
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profile as a result of dynamical friction (e.g., Dobke et al. 
|2007| |Auger||2008| ). We therefore directly probed the slope 
of the mass prohle of quasar host galaxies to test whether 
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quasars do in fact tend to form in steep-profile dynamically 
perturbed galaxies or in shallower-profile (isothermal) more 
relaxed galaxies. 

We compared the mass density profile of PSQ host 
galaxies with two comparison samples. We first compared 
them to a sample of galaxies that did not harbor an active 
central black hole and a second that did but that did not 
show any evidence for recent merging activity. For the first 
comparison sample we considered 22 non-quasar early type 


galaxies from the Sloan Lens ACS Survey (SLAGS; Bolton 


et al. 20061 for which weak and strong lensing measurements 
had been made ( |Gavazzi et al. ] |2007| hereafter G07). The 
second was a sample of quasars from the Co-evolution of 


spheroids and black holes HST program (Treu et al. 2007 


hereafter T07). These moderately luminous Seyfert galax¬ 
ies represented a typical set of quasars whose host galaxies 
did not have strong star formation but did contain an ac¬ 
tive galactic nuclei ( [Barger et al.|[^14[ ). By comparing the 
halo of the PSQs to the SLACS and T07 halos we examined 
any differences that may indicate a connection between the 
mergers that triggered a quasar and the mass profile of the 
quasar host galaxies. 

Tracing the differential mass components of a galaxy 
has now become commonplace amongst large scale sky sur¬ 
veys. Exploiting a galaxy’s gravitational influence on the 
local space-time allows a non-parametric and assumption- 
free analysis of the projected density profile. Weak gravita¬ 
tional lensing of galaxies by galaxies or ‘galaxy-galaxy leas¬ 
ing’, describes the effect that a foreground galactic potential 
has on the observed shape of a background, ‘source galaxy’. 
By measuring the shape of background galaxies and corre¬ 
lating them to foreground galaxies it is possible to extract 
the underlying projected density profile. The first detection 
of galaxy-galaxy lensing was by made by [Brainerd et al.j 
( |1996| ). This pioneering analysis studied the gravitational 
lensing signal around galaxies and attempted to fit a sin¬ 
gle component mass model to the data, however this model 
assumed all the mass was contained within the galaxy halo 


(Schneider fc Rix|1997 Hudson et al.|[l9M Hoekstra et al. 


2004). Subsequent studies found that either analysing iso¬ 
lated galaxies ( [Hoekstra et al.||2005| ), or applying more so¬ 
phisticated models that included second halos were required 


(Mandelbaum et al.|2006' van Uitert et al.|2011 Leauthaud 


et al. 20121. One of the most recent study was from the 
CFHTLenS survey ( Heymans et al.|2012 l. Studying the re¬ 
lationship between the dark matter and stellar halo around a 
galaxy ( Velander et al.|2014 1, this ground based survey used 
galaxies at 0.2 < 2 : < 0.4 and applied a multi-component 
model to dissect the contributions from stars, dark mat¬ 
ter and surrounding satellite halos. Although complimen¬ 
tary to the study here, the advantage of 155 square degree 
survey was that the analysis could extend out to 10 Mpc, 
whereas here we analysed the more central parts, in the 
range 10 < r < 1000 kpc, where quasar feedback has most 
impact. 

Studies specihcally probing the host galaxy of a quasar 
are few. Most recently, the profile of galaxies hosting Active 
Galactic Nuclei (AGN) has been mapped using the COS¬ 
MOS quasar catalogue ( Leauthaud et aL][2014 l. At a mean 
redshift of 2 = 0.7, they used 382 moderately luminous X- 
ray AGN and found that quasars tended to lie in medium¬ 
sized halos with Mhaio ~ lO^^M©, with a halo occupation 


distribution not too dissimilar from normal galaxies that do 
not host quasars. At this redshift the majority of the AGN 
are fueled by the galactic bar and therefore exhibit a similar 
profile to that of normal galaxies, which is indeed supported 
by their measurements. Prior to this, [Mandelbaum et al.[ 
(20091 studied the dark matter halo of radio loud quasars in 


the Sloan Digital Sky Survey and compared them to a sam¬ 
ple of control galaxies Hnding that these AGN laid in larger 
dark matter halos than the galaxies. In this study we build 
upon and extend this body of work by directly probing the 
dark matter halo of galaxies that are in a specific period of 
their life, where discrepancies in the dark matter halo may 
be most evident. 

In the following, we adopt a flat cosmology with H = 
100h“^km s“^Mpc“^, Dm = 0.3 and Ha = 0.7. 


2 DATA & METHOD 

The HST snapshot survey of 29 post-starburst quasars (Pro¬ 
posal ID 10588; Gales et al.pOll 20131 was well suited to 
our study. The sample was imaged using the HST F606W 
filter on the Advanced Camera for Survey (ACS) at a mean 
redshift of 2 = 0.32 . We compared these galaxies to two 
control samples built from archived HST/ACS programs. 
The hrst was composed of massive elliptical galaxies from 
the Sloan Lens ACS Survey (Proposal ID 10494 and 10886), 
for which 22 strongly lensed galaxies were imaged through 
the F814W filter of ACS ( [Bolton et al.|200^ G07) at a mean 
redshift of 2 = 0.23. These galaxies had been shown to be 
representative of the general populations of elliptical galax¬ 
ies ( Barnabe et al.|2011| . The second contained 17 Seyferts 
galaxies that hosted moderately luminous AGNs at a mean 


redshift of 2 = 0.36 (Proposal ID 10216 Treu et al. 2007). 
The host galaxies of these quasars showed no obvious sign 
of any merger with no signihcant active star formation. For 
detailed selection criteria we point the reader to the cited 
references. 

In all cases we first corrected the images for charge 
transfer inefficiency ( Massey|2010 Massey et al.|20T4 l, and 
we calibrated and drizzled each image individually for PSF 
estimation and then we combined them for our final shape 
measurement. We then used the RRG pipeline from the 


COSMOS weak lensing survey (Rhodes et al. 2000 Leau¬ 


thaud et al. 2007) and developed by Harvey et al. (2015), 
to measure the shapes (ellipticities) of the faint galaxies in 
the field of view. Although we did not have photometric 
or spectroscopic redshifts of the background galaxies, the 
quasars were at sufficiently low redshift (2 < 0.5) that we 
could assume that all measured galaxies were in fact behind 


and at a redshift of 2 = 1.0 (Leauthaud et al. 2007). We 


initially made the same magnitude cuts as G07, removing 
galaxies that were outside the interval 20 < M < 26.5 how¬ 
ever we found that we obtained significant contamination 
from bright sources near the lens and therefore cut between 
21 < M < 26.5 in all three samples. We also cut in size 
removing any galaxies that are less than 3 pixels in size and 
a signal-to-noise of less than 4.4. We also removed pairs of 
galaxies that were within 0.5” of one another. 

Having measured the ellipticities of the background 
galaxies, we azimuthally averaged the tangential ellipticity 
in radial bins around the foreground objects in the three 
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above samples and we converted this to the projected sur¬ 
face density in a similar way to G07 via 


AE = 






—2 

t,j 2^i=l 


E Nlens 
J = 1 


y-2 


e :: 


( 1 ) 


where AE is the average excess surface density in a given 
radial bin, 

fp n 

Scrit = (2) 


47rG DisDoi ’ 

et is the tangential ellipticity around the galaxy we are prob¬ 
ing, and (Je is the error in the ellipticity, and subscripts re¬ 
fer to the relative angular diameter distances between the 
observer(o), the source(s) and the lens(l). To calculate the 
value of Ecrit for each lens we used the reported redshift and 
assume that the source is at a redshift of Zsource = 1-0. The 
hrst summation in equation Q is over N galaxies that lie in 
the radial bin, and the second summation is over the sample 
of galaxy lenses, Aliens- 


3 RESULTS 

The top three panels of Figure show our results for the 
different samples of galaxies, where errors shown are the Icr 
error on the mean. Given the small sample we jack-knifed 
each lens to ensure that no one object was biasing our re¬ 
sults. We found the errors derived from this test were con¬ 
sistent with those shown here. We fitted a Navarro, Frenk 


the large-scale dark matter halo (|Bartelmann 

|l996|l and a 

de Vaucouleurs profile (de Vaucouleurs 1948 

Maoz & Rix 


19931 to model the more central part of the total mass pro¬ 


file. Note that in doing so, we did not attempt to use any 
prior knowledge from the light distribution other than fix¬ 
ing the effective radius. This would require assumptions on 
the stellar populations of the galaxies, which were all very 
different in the three samples we use. We restricted our¬ 
selves in the present work to the total radial mass profile 
using galaxy-galaxy weak lensing. We fixed the effective ra¬ 
dius of the de Vaucouleurs in each case using the reported 
results from the three main papers and hence only fitted 
the normalization of the de Vaucouleurs, and the normali¬ 
sation and scale radius of the NFW halo. In each panel of 
Figure we show the best fitting model and in grey the 
two contributing components. Tablegives the best htting 
masses derived from the fits in units of Mq, the de¬ 

rived NFW concentration and the chi-square of the best fit 
over the number of degrees of freedom. As a cross refer¬ 
ence, we found that the NFW masses for the SLAGS sample 
were consistent with the NFW masses in G07 analysis of 

AfNFW = 14^5 X 

The result derived from the two quasar samples im¬ 
plied they were compatible and drawn from a similar pop¬ 
ulation. To test this conclusion we compared their velocity 
dispersions. Using the known relation between black hole 


mass and velocity dispersion from Tremaine et al. (20021 


we found that the PSQ sample had a mean velocity disper¬ 
sion of = 208 ± 41 km/s and the T07 sample had 

a velocity dispersion of = 228 ± 36 km/s. Given the 

small sample size, we carried out a KS test to decide whether 
the two samples could have been drawn from two different 



Radius [Mpc/h] 

Figure 1. Weak lensing signal for different types of galaxies. 
From top to bottom we show the excess projected average surface 
density for: a) the 28 host galaxies of PSQs, b) 22 massive early- 
type galaxies from SLAGS with no active supermassive BHs at 
their centre, c) the 17 Seyfert galaxies of T07 and d) the combined 
quasar sample T07+PSQ. Note that for the SLAGS sample we 
perform our own weak lensing analysis. Our measurements are 
consistent with the ones from G07. The shaded region denotes the 
regime that is too close to the lens to gain robust weak lensing 
constraints, and hence any fit is a direct extrapolation. 


Sample 

AL/auc 

Mnfw 

C200 

X^/dof 

PSQ 

1.7± 1.2 

4.0 ± 3.0 

3.6 ± 2.8 

4.90/5 

SLAGS 

< 0.84 

13.0 ±3.6 

5.4 ± 1.5 

4.35/5 

T07 

0.72 ± 0.8 

4.4 ± 3.2 

4.1 ± 3.0 

2.51/5 

PSQ + T07 

0.18 ±0.57 

2.8 ± 1.9 

4.9 ± 3.5 

5.75/5. 


Table 1. The derived masses and concentration parameters from 
the best fitted models to the data. The first column stipulates 
the sample of galaxies, the second column gives the mass within 
the effective half light radius of the best fitting de Vaucouleurs 
profile in Mq, the third column gives the best fitting 

NFW mass also in Mq and the fourth column gives the 

derived concentration with error. The final column gives goodness 
of fit to data for the reported model (chi-square / degrees of 
freedom). 
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populations and we found no significant evidence that they 
did {P — 0.35). Given their similarity we combined the two 
sets of quasar host galaxies to test whether the resulting 
increased sample differed from the SLAGS sample. 

The bottom panel of Figurej^gives the combined profile 
of galaxies that contain an active supermassive black hole at 
its centre. We found that when combined there was no sig¬ 
nificant difference in the profile shape between these galax¬ 
ies that did harbour an active central black hole and the 
profile of a galaxy that did not (SLAGS). In addition, this 
profile did not significantly depart from that of an isother¬ 
mal sphere. This result was consistent with |Leauthaud et al.| 
120141 in which they found no difference in profile for galax¬ 
ies containing moderately luminous quasars. However we did 
find a significantly different NFW mass, with the quasar 
sample lying in much smaller dark matter halos than the 
SLAGS galaxies. We compared the velocity dispersions of 
the two quasar samples with the SLAGS sample and found 


that the SLAGS sample with (aj 


=) = 244 ± 44 km/s 


was consistent with the T07 sample (P=0.59) however, the 
distribution of velocity dispersion did differ from the PSQ 
by ~ 2(7 (P=0.05). Our conclusion is that there is a differ¬ 
ence between the PSQ host galaxies and the SLAGS host 
galaxies, however this cannot be attributed to recent merg¬ 
ing activity and is most likely due to the environment and 
formation history. 


3.1 Systematic Tests 

We verified our weak leasing measurements by carrying out 
two standard validation checks. We first checked the ex¬ 
cess surface density derived from the cross-shear, i.e. the 
stretch along the axis 45° to the tangential shear. Under 
the assumption there are no other external factors such as 
lens light impacting the shape measurement or poor PSF 
removal, the cross-shear should average out to zero, and 
hence the surface density calculated should also be zero. Our 
second test measured the tangential shear around random 
points in the field, which should also be zero since they are 
not correlated in any way. We found in both cases that all of 
the samples exhibited no significant excess signal. We also 
carried out a similar test but with the two components of 
ellipticity: el and e2, where in the absence of any additive 
bias the mean of these will be equal to zero. We check all 
three samples and we found no evidence for any significant 
bias. 

In addition to the cross-shear we also carried out checks 
to test for bias due to lens light, source galaxy redshift as¬ 
sumptions, noise bias and source galaxy blending. We first 
checked that the presence of lens light was not a significant 
issue with our lensing estimation. To do this we measured 
the median lens light profile for each sample and normalized 
with respect to the intensity at a radius of 1 arc-second 
(where /(i?=l")psQ = 0.046 e“/second, I{R—l")'Tm = 
0.028e“/second, /(i?=l”)sLACs = 0.17e“/second). The top 
panel of Figure shows that the two quasar samples, PSQ 
and T07 both dropped to the sky background by ~ 4” or 
~ 15kpc and hence did not impact our shape measurement. 
The SLAGS lenses were slightly larger, extending out to 
~ 8”, which could have impacted our source selection and 
shape measurement. To ensure that this was not the case we 
carried out two further checks. We first measured the signal- 



Radius [arcseconds] 


Figure 2. Systematic tests on the different samples of galaxies. 
We show in the top panel, ranging from 0.01 to 10 arc-seconds, 
the light profiles of the three different samples. The two quasar 
samples are comparable, and the contribution from the quasars 
can be clearly seen. The middle panel shows the median signal-to- 
noise for each source galaxy used in the analysis, and the bottom 
panel shows the source galaxy density as a function of distance 
from the centre of the galaxy. The middle and bottom panels span 
the entire field of view from 7 - 100 arc-seconds. 


to-noise ratio (S/N) over the field for each sample, checking 
to see if sources were detected within the lens (and hence 
have unreliable shapes). Such an effect would have mani¬ 
fested itself as a mean increase in source galaxy S/N. We 
also examined the number density of galaxies over the field 
of view, checking to see if the source selection had picked up 
substructure within the lens that could mimic source galax¬ 
ies. The bottom two panels of Figure show the results. 
We found that the S/N was constant across the field for all 
samples suggesting a robust source selection, and we also 
reported a constant number density of galaxies for the two 
quasar samples, also suggesting we did not select substruc¬ 
ture within the lens. As expected the number density for 
the SLAGS lenses drop in the inner most bin, this is due to 
an increase in noise due to foreground larger SLAGS lenses. 
In addition to this we found that the innermost bin of the 
cross shear for the SLAGS data was consistent with zero, fur¬ 
ther supporting the robustness of the source selection and 
shape measurement. The S/N also shows that the galaxies 
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selected had a mean S/N> 15 for all samples, where noise 
bias is negligible ( Amara fc Refregier|2008 1. 

The bottom two panels of Figure show the inter¬ 
sample variation in the number density and S/N. This sug¬ 
gested that the mean source redshift also varied between 
samples. Since we assumed that all sources were at a red- 
shift 2 = 1.0 we tested the impact of this assumption on our 
results. We re-ran the analysis with a variety of different as¬ 
sumed source redshifts and studied the effect on our results. 
We note here that changing the source redshift would not 
affect the profile shape of the measured galaxies but only 
the mass estimates as the source redshift only affects Ecrit 
in equation (§. However, despite this, we found that if we 
varied the source redshift from 2 = 0.8 : 1.5, the change in 
derived masses were negligible compared to the statistical 
error bars. We therefore ignored this source of error. We fi¬ 
nally checked to see if results were robust to source galaxy 
blending. In the method we removed all pairs of galaxies that 
were separated by less than 0.5”, to check that our results 
did not depend on this we removed all galaxies that were 
separated by less than 1” and we repeated the experiment 
with an angular separation of 1”, with the same results. 
We found no significant change in the profile of the galaxies 
and the derived masses, and were therefore confident in the 
source selection. 


3.2 Future surveys 


Our study of 29 PSQs showed that the dark matter halo in 
these galaxies did not depart from isothermality. Given the 
accuracy of our measurements we tested how sensitive our 
measurements were to a change in profile, and how much it 
would need to change in order for us to detect a discrepancy. 
Additionally, with the advent of all sky space surveys we in¬ 
vestigated how sensitive future measurements are expected 
to be. To do this we fitted a very simple power law profile to 
the data, AE oc r“ and measured the index, a and its asso¬ 
ciated error. We then artificially steepened the data, raising 
the index and refitted the profile, measuring the index with 
an associated error. We then compared this to an isothermal 
profile ( Q = — 1) and determined how significant the depar¬ 
ture from isothermality we had measured. As the data be¬ 
came steeper we were able to better detect a departure from 
isothermality. Figurej^shows the relation between the mea¬ 
sured gradient of the profile and the measured significance 
of a departure from isothermality. We repeated this exercise 
using different fake sample sizes, decreasing the error bars 
by a factor -^/Aquasars. We found that with future surveys we 
will be able to detect a 10% departure from isothermality at 
the la confidence level using ~ 1000 quasars. Such a sample 


size will be possible with surveys such as Euclid (Laureijs 
et al.|20TT I in the near future. 


4 CONCLUSIONS 

We measured the total mass profile based solely on post 
starburst quasars (PSQ) host galaxies that showed signs of 
recent merger activity. Our measurement, based on weak 
gravitational lensing of background galaxies, is compared 
with that of the SLAGS sample of massive elliptical galax¬ 
ies and with a sample of luminous Seyfert galaxies. We found 



Figure 3. The sensitivity of our data to a change in profile. We 
artificially steepened the data to create a fake ‘steeper’ profile, 
and fit a power law profile, measuring the index a, and the as¬ 
sociated error. We found that for our sample size of 29 quasars, 
we would measure a la deviation from isothermality if the profile 
had exhibited an index of a ~ —1.2 (where Oiso = —1). Future 
surveys will observe orders of magnitude more quasars resulting 
in more sensitive data. The dotted black line shows a detected Itr 
departure from isothermality. 


that in all three cases the profiles did not depart significantly 
from isothermallity and concluded that within the error bars 
shown, mergers do not affect the shape of the dark matter 
halo profile. Given the commonality of the two quasar sam¬ 
ple in NFW mass, de Vaucouleur mass, and central black 
hole mass, we found that AGN activity did not affect the 
shape of the dark matter profile, and therefore in these two 
samples of galaxies baryonic feedback was not strong enough 
to affect the potential via gravitational or non-gravitational 
forces. 

This study has presented the first total mass and mass 
profile measurement of non-obscured, low redshift quasars 
that show recent signs of merging. Given the consistency 
of the results and the small sample size we tested to see 
how sensitive the data actually was to a change in profile 
shape. We found with the current data we would be able 
to detect a ~ 20% departure from isothermality at the Icr 
level, however with future surveys we will have access to 
many more quasars, reducing this sensitivity to a > 10% 
change in profile shape. However, these estimates are purely 
statistical and the robustness of these results would need to 
be tested in the event of more accurate data. 
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